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Abstract: We show that in supersymmetric grand unified theories new effective D=4 
and D=5 operators for proton decay are induced by soft SUSY-breaking terms, when 
heavy GUT gauge bosons are integrated out, in addition to the standard D=6 ones. As a 
result, the proton lifetime in gauge mediated channels can be enhanced or even suppressed 
depending on the size of the heavy Higgses soft terms. 
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1. Introduction 

One of the immediate consequences of Grand-Unified Theories (GUT) is baryon number 
non-conservation that can lead to proton decay jj], H. The heavy gauge bosons mediate 
the effective baryon-number violating four-fermion operators 

O gaU ge ~ g 2 G M G 2 (qu c Jd c ), g 2 G M G 2 (qu c qe c ), (1.1) 

where Mq and g G are the Grand Unification scale and gauge coupling constant, q = 
(u,d)L, I = (y,e)L ar e the left-handed quarks and leptons (weak isodoublets) and u c , d c , 
e c are charge-conjugated fields of the right handed ones: ur, c£r, ej? (weak isosinglets) . In 
addition, other four-fermion operators of different structure are mediated by heavy colored 
triplet Higgses with mass Mh ~ M G : 

Ohiggs - g 2 Y M H 2 (q qql), g 2 Y M H 2 (u c u c d c e c ) , (1.2) 

The latter operators are typically weaker because of smallness of the Yukawa couplings gy 



but in some models they can be dominant over the gauge mediated operators (1.1 



The set (1.1), (lO) represents all possible D=6 baryon number violating operators, 



independently of the details of grand unification pi. The two kinds of operators have 



different chirality structures (LRLR for (1.1) and LLLL or RRRR for (|1.2|)) and they 



could in principle be distinguished via the polarization of final states y]. Both kinds are 
suppressed by two powers of Mq. 

In non supersymmetric models, when Mq is below 10 15 GeV, processes following 
from (|1.1| ) are already ruled out. Supersymmetry, in addition to making the unification 
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more natural, raises its scale, setting the magnitude of these processes within the reach of 
near future experimental facilities. 

Supersymmetry raises the unification scale to 10 16 GeV and makes these operators 
hardly observable. On the other hand it introduces additional D=5 baryon number violat- 
ing operators mediated by heavy colored higgsinos ||, as 

Ohiggsino ~ gyM^iqqql), g\ M H l \u c u c d c e c ) , (1.3) 

where tilded fields represent scalar superpartners. These are suppressed by a single power 
of the GUT scale and after dressing by gaugino exchange they give rise to operators of 
the form (1.2) with a cutoff scale ~ (M#ms) 1//2 where is the SUSY-breaking scale ||. 



They thus become generically dominant and on the verge of being in conflict with current 
experimental limits on these specific decay modes (p — > Kv etc.) [Q]. However, their 
magnitude is very model dependent: essentially they exclude minimal versions GUTs and 
cause problems for models unless fine tuning is arranged. Several mechanisms have been 
devised to suppress them by playing with the structure of the heavy sector of the theories. 1 

Let us remark that gauge coupling unification does not strictly require supersymmetry 
of the theory. For instance the presence of fermionic partners of the gauge and higgs at 
TeV scale can adjust the running of the gauge coupling constants so that they unify at one 
point. Though scalars are not crucial for unification they are predicted by low scale SUSY; 
however, finding at LHC a SUSY-like spectrum would not mean that supersymmetry is 
discovered. Indeed, it would be extremely difficult to verify that the lagrangian has a 
supersymmetric structure, i.e. that the different coupling constants are related, like the 
quark-squark-gluino coupling constants that should be exactly the same as the strong 
gauge coupling constant. 

One can imagine a fake-SUSY theory where only the sparticle spectrum is supersym- 
metric (i.e. every particle has its "superpartner" ) while the the lagrangian is not. What 
would happen in such a theory? We argue that, even if the gauge coupling unification is 
achieved as perfectly as in a truly supersymmetric theory, it would lead to disastrous pro- 
ton decay rate. The reason is the following: once such a theory contains scalars partners 
of quarks and leptons (q, I) it generically contains D=4 operators of the form 

O quartic ~ {q*u c l*d c ), (q*u c q*e c ), (1.4) 

which in a GUT context can not be excluded by any symmetry reason. Notice that even if 
they are not present in the bare lagrangian they emerge radiatively by loops of GUT gauge 
bosons. The dressing by gauginos transforms these D=4 into D=6 ones on the form ( |1 . 1| ) 
that directly cause the proton to decay at a dangerous high rate, being suppressed only by 
two powers of the fake superpartners mass scale that is of order TeV. 



1 There are several ideas how dimension-5 operators can be suppressed. In particular this can be due to 
special arrangements in the heavy higgs sector || , because of symmetry properties of the Yukawa sector Q . 
In SO(10) models, LLLL operators can be naturally suppressed by the choice of the SO(10) breaking VEVs 
while the less dangerous RRRR ones are left allowed. With further model building also these latter can 
be eliminated Finally, in supersymmetry there are also D=4 B and L violating operators that can be 
forbidden by exact R-parity 0. 
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Complete supersymmetry instead provides an automatic protection from these D=4 
operators: they in fact correspond to the D-terms relative to the broken gauge generators, 
and they have to vanish if SUSY is unbroken. What happens is that the existing D-term 
involving light fields g 2 \4>*(j)\ 2 is cancelled by two other diagrams: one with the exchange of 
the (broken) gauge field and one with the exchange of the heavy longitudinal part of the 
higgs field that breaks the gauge group. For this cancellation to hold it is crucial that the 
coupling g is exactly the same in the three graphs, i.e. that SUSY is exact. The shortest 
proof of this fact can be given in the superfield formalism, where the only supersymmetric 
D-term involving four light superfields is [<^<l><I>t<3?]£). This operator does not contain a 
four scalar contact interaction, that therefore has to vanish. 

However, since supersymmetry has to be broken, one expects that this protection 
mechanism works only partially and that the susy-breaking terms will turn on such op- 
erators. As a result they may significantly affect the proton decay. Indeed, it was noted 



in [10] that SUSY-breaking induces the D=4 scalar operators (1.4). However, surprisingly 
enough a complete analysis of the soft-susy breaking effects on proton decay has not been 
performed. 2 

In this work we study the effect of the soft terms on the low energy effective theory 
produced after the heavy gauge superfields are integrated out at the GUT scale. We 
show that the D=6 operators are always accompanied by new operators of D=5 and D=4, 
turned on by the presence of the soft terms. Next, we compute the renormalization of 
these operators from the GUT scale to the SUSY breaking scale; we adopt the techniques 
illustrated in |0| that simplify considerably the task. We then dress the new D=5 and D=4 
operators at the SUSY breaking scale, transforming them in the form (|1.1[) and estimate 
when they can be relevant. As an example, we discuss the SUSY SU(5) model and show 
that their contribution can be important and could bring the proton decay rate in specific 
channels to be experimentally accessible. 

2. Gauge mediated effective operators in softly broken SUSY GUT 

The gauge mediated effective operators are efficiently described in the superfield formal- 
ism with soft breaking terms inserted as spurions. If we arrange the chiral superfields of 
irreducible representations in a column vector $ = {<&/}, the full lagrangian is, in compact 
notation: 

C = Jd A 6 &Xe 29V <5> + Jd 2 9 W($) + W a W a Y +h.c. (2.1) 

Here the gaugino masses enter via Y = (l+rrigO 2 ), and the soft D-terms may be parametrized 
by the matrix X = (1 + T8 2 + TW 2 + Z9 2 6 2 ) via the matrices {T u } of order m s and {Z LJ } 
of order m|. One can however perform a field redefinition to set Tjj = 0. For simplicity 
we will also consider only universal soft terms, taking Zu diagonal. Therefore we have: 

X I j = X I 5 I j = (l-mje 2 e 2 )5 I j. (2.2) 



2 In it is described a classification of all the D=4, 5, 6 operators relevant for proton decay, while in 
| the effect of soft terms in a SUGRA scenario was studied, but only for the analytic D=5, 4 operators. 
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The superpotential W(3>) includes the soft F-terms via spurion fields and may be parame- 
trized similarly (see e.g. fl4.f|) ) but its explicit form is not directly relevant for this section. 

At the scale of gauge symmetry breaking one can decompose in (&h, </>%), 
respectively heavy superfields, light goldstone superfields, and light non-goldstone super- 
fields. The decoupling of heavy superfields (e.g. colored higgses) leads to dimension-6 
and analytic dimension-5 effective operators that may violate baryon number. The decou- 
pling of heavy gauge fields and goldstones in turn leads to the D-term effective operators 
of dimension 6 All these operators are affected by the soft susy breaking terms. 

To find the effect on the gauge mediated dimension-6 operators, it is convenient to 
adopt the so called super-unitary gauge jl4]l , where the goldstone superfields <3?a are gauged 
away inside the broken massive gauge superfields, denoted as V A . To integrate out Va one 
expands the gauge exponential in (^) to the quadratic order 

+ 2J A V A + KabVaVb + ■■■ 
J a = g&XT A $, 

K AB = 2g 2 &XT A T B <S> (2.3) 

and then one notes that in the unitary gauge 

J A = 9 4 T A& Xi K AB = 2g 2 (Qh^TaTb ($ h ) X h . 

As expected a VEV of the heavy fields give a squared-mass matrix Kab to the broken gauge 
fields. 3 In a suitable basis of broken generators this matrix is diagonal, Kab = Ka$ab- 
We can note already at this stage that the heavy gauge boson mass matrix contains SUSY- 
breaking factors, such as the X's. 

The result after integrating out the broken gauge fields Va is then: 

J^OJaKiIJb = 

-hwwA^^i^) (^) • <") 

where summation on all indices is understood. Note that in the integration we have ignored 
sub-leading terms like the gauge kinetic term and the gaugino masses for the V A gauge 
fields. In fact at they both give subleading effects in fl2.4fl. 

Considering that ($h) ~ Mq, we recognize in the first factor the standard coupling 
constant of the dimension-6 operators ~ 1/Mq. The supersymmetry breaking however has 
propagated in this operator, and indeed the second factor involves the soft SUSY breaking 
D-terms X that are carried along in the decoupling process. Moreover, due to the soft 
SUSY breaking in the superpotential, also has in general a non-vanishing F-term, 

(&h) = ^ff (1 + IhG 2 ), that induces an other supersymmetry breaking in the effective gauge 
bosons mass. 

In the presence of non-universal soft terms Ja has an additional piece {^HyTAXH^j, that gives new 
soft masses to (f>p Also Kab is modified in a similar fashion, see [[[3]]. However the present analysis is not 
affected substantially. 



£ (2) = d 4 e 
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Figure 1: New operators generated by soft SUSY breaking terms in the heavy gauge exchange. 




Figure 2: four-fermion operators after dressing via gluino exchange. 
The overall result with SUSY-breaking terms can be conveniently rewritten as: 

£ I A 6 (l + £9 2 + ^9 2 + oj^ 2 ^ U\T A< k\ (^)Ta^ 3 ) d 4 9 • (2.5) 
ij,A J 

where A6 = g 2 /M\ is the super symmetric four-fermion coupling, the heavy gauge-bosons 
masses are given by M 2 A = J2h^9 2 ( v h^aTavh), and the SUSY-breaking coefficients are: 

£ = -f H , cuij = -m 2 - m 2 +m 2 H + \f H \ 2 • (2.6) 

For simplicity in this last expression we have assumed the breaking by a single VEV. 4 

In terms of field components the effective operator fl2.5|) contains the three operators 
shown in figure [lj the standard dimension-6 four-fermion operator ijj*ipijj*ip with coupling 
~ X/Mq\ then a new dimension-5 operator of the form A*il)A*il) + h.c. with coupling 
~ tus/Mq coming from the terms with 6 2 and 9 2 , and finally a new dimension-4 operator 
of the form AA*AA*, with coupling ~ m 2 s /M^. 

The new dimension-5 and dimension-4 operators can be dressed by gaugino exchange 
at the SUSY-breaking scale (see figure |2|) and transformed in effective dimension-6 four- 
fermion operators, as it happens for dimension-5 analytic operators. Each dressing loop 
brings a factor ~ 1/ms, so that the effective strength of all these operators is the same, 
1/Mq. The actual relative strength will depend on the coupling constants involved in the 
dressing and on the ratio of the effective soft breaking parameters £, u to the gaugino 
and/or sfermion masses. 

4 With more VEVs, in the first formula fn should be replaced by its "average" 
(Ejf^l(if)/ fl )/(EH^l(H))i where M\i H -, = 2q 2 v' h TaTavh ■ Similarly in the second formula for 
m 2 H and \/h\ 2 - 
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3. Running and dressing 



To calculate the effect of the three operators in (2.5) one has to run them from the decou- 
pling (GUT) scale down to the SUSY breaking scale and dress them to get the dimension-6 
effective operators. The running below the SUSY scale is non-supersymmetric and was 
analyzed in [16]. 

Renormalization mixes the supersymmetric and the non-supersymmetric effective op- 
erators via the soft susy breaking parameters of the theory, mainly the gaugino masses. The 
detailed computation is rather complicated due to the large number of diagrams involved. 

Instead of attacking the problem by brute force, we employ the elegant techniques 
devised in to analyze the soft terms renormalization. Starting from the anomalous 
dimensions of the supersymmetric operators, we find the renormalization in the softly 
broken theory by promoting the couplings to full superfields built with the soft terms. 

We start from the definition of the renormalization of the supersymmetric coupling 
in (l2l 



A 6 



a 3 , a 2 , a.\. 



n 

i=l,2,3 



or 

a; 



(0 
6« 



(3.1) 



where &w is the beta-function coefficient for each gauge group and 
ing supersymmetric contribution to the anomalous dimension of A 
dimensions were calculated in [[l?]] 



(3) 
76 



4/3, 7 < 2) = -3/2, 7e 



(i) 



7g^ is the correspond- 
ed These anomalous 
-23/30. 



The key step, to renormalize the full operator Xq(\ + !^Q 2 + !^6 2 +u)Q 2 Q 2 ) in the presence 



of soft terms, is to promote each gauge coupling a to a 

xj(i + ( B e 2 + ^ B + to B e 2 e 2 ) _ 
a 6 (i + ^ 2 + ^e 2 + uj9 2 9 2 ) 



a(l + m g 9 2 + mW 2 + 2\m g \ 2 6 2 6 2 ): 



Z e (a 3 ,a2,ai] 



(3.2) 



Expanding then Z§ in grassmann variables we find how the operators of different dimension 
mix under renormalization: 



47T- 



df 



/A 6 \ 

A 6 £ f 

\X 6 UJ ) 



( 1 



76« 



m g 




1 






1 

17ln 



\ 



/A 6 \ 

A 6 e f 

\X G UJ ) 



(3.3) 



\2|m g | 

where t = ln(/i 2 ) and a summation on the different gauge groups is implicit in the r.h.s.. The 
gaugino masses m g and gauge couplings a follow the equations m g /m g = a/a = bet /An. 



J 



The equations (|3.3| ) are solved in terms of the evolution of the gauge coupling constants 
a from the GUT to the SUSY scale by using the auxiliary functions 

a{S) 



R 



a(G) 



R 



76 (R 



i 



Ri 



76 ( r>2 



i 



In the one-loop approximation only the renormalization due to gauge loops needs to be taken into 
account, and in the operators involving the first generation the contribution of large top Yukawa is sup- 
pressed by mixings. In higher loop orders one should also include other effects, for example the threshold 
corrections due to insertions of more than one A6. 
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The result is: 



\ 6 (S) = IT^\ 6 (G) 
^S) = aG)+Ri m~ g {G) 

uj{S) = uj{G) + 2i2i Refc(G) m|(G)] + (i? 2 + flf) mJ(G) . 

For example, the largest effect comes from SU(3) color, for which 7^ = —4/3, = —3, 
a 3 (M z ) = 0.119 and a 3 (G) = 1/23: 

A 6 (S)~1.55A 6 (G), (3.4) 
£(S)^(G) + 0.74 m^(G), (3.5) 
w(5) ~ w(G) + 1.48£(G) m H (G) + 3.26m? 3 (G) , (3.6) 

where for simplicity we assumed £ and real. 

The effective strength of the dimension 6, 5 and 4 operators can be compared after 
dressing with exchange of some gaugino. As shown in figure ^, it is clear that the chiral 
structure of the D=5 operators requires a Majorana mass to perform a chirality flip, and 
for low momentum processes as proton decay this is true also for the D=4 operators. The 
D=5 operators can be dressed by gluino exchange; the D=4 operators on the other hand 
can only involve one gluino exchange while the other loop is necessarily formed via W-ino 
(or B-ino exchange). 

We can estimate the strength of the two new operators by defining the corresponding 
effective four-fermion couplings at the SUSY scale: 



A 5 = A 6 2£ — L(m i ,mj/,m §3 ) 



A 4 = A 6 oj ^ L(mx,my,mg 3 ) L(mx',my>,mg 2 ) 

where the sfermion masses entering the dressing loop(s) and a 2 , «3 the gauge 

coupling constant involved. All quantities are evaluated at the SUSY scale. L is the loop 
integral: 

2 2 i , 2 2 i m o , 2 2 i m \ 

mfm.2 log —j + m^m^ log — f + mf m§ log — f 



[nil — in-,) (in- — nijj (pin — ni'i 



L(m 1 ,m 2 ,m 3 ) = m 3 - , ., Wi ., ,- 2w ■> 



n ~~ '"-2; i"n ~~ '"-3; 
i a* - 1 - log 

= z 7) for mi = ?7i2 = m (3.7) 

777,3 (ZKr — l) 

The loop integral is plotted in figure ||, where all the masses are measured in TeV. From 
there we see that L may be of order ITeV -1 (with a maximum of L ~ 5-6 TeV -1 ) for 
small sfermion masses ~ 100 GeV. On the other hand the gaugino mass may be raised up 
to 1-2 TeV before starting to suppress the loop. 
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Figure 3: Loop function L in TeV , with m§ ranging from 0.1 (upper) to 1 TeV (lower) in steps 
of 100 GeV. 



The present limits [27] allow squark masses as low as 100 GeV when the gluino mass is 
> 500 GeV. In this parameter region, we find that L is almost maximal, L ~ 5, so for the 
numerical estimates we will stick to this choice. A different choice can be easily considered 
by extracting the relevant loop factor from eq. ( |3.7|) or directly from figure ||. 

With this choice the strengths of the new operators relative to the D=6 one are: 

A5/A6 ^ — 10 ■ 



4vr TeV 10 TeV 

A4/A6 "4^4^ 25 TeV"(3b^vF- (3 ' 8) 

As a result, the effect of D=5 and D=4 operators may be comparable (or larger) than that 
of the standard D=6 operators. However for this to happen the effective susy-breaking 
terms £ and ui should be larger than the soft masses. One needs for example £ ~ 10 TeV, 
a factor of 20 or 100 larger than the gaugino or sfermion masses. 

One should also ask whether these large values might be generated in the evolution of 
13 orders of magnitude from the GUT down to the SUSY scale, by mixing with other soft 
parameters, namely the gaugino masses. However from the running ( |3.6|) we see that in 
the regime of £, uj 3> rrig 3 the gaugino gives a small contribution that does not modify the 
estimate (|3.8|). 

Is it then plausible for £ or uj at GUT scale to be so larger than other soft susy breaking 
parameters in the theory? We argue that this is possible without spoiling the framework 
of low energy supersymmetry. The reason is as follows: from the expression of £ and u>, 
eq. ( |2.6D , we see that they are induced, in addition to the soft masses of fermion fields, 
by m 2 H and fn, the soft SUSY-breaking parameters in the heavy higgs sector. One can 
not play much with the soft masses of the heavy fields rrijj, since these are constrained 
because they usually mix with the MSSM higgses soft masses in the renormalization from 
the Planck to the GUT scale. On the other hand the F-terms fjj are less constrained, since 
they do not directly enter in the running and one should not assume them to be small. 
This can be seen in the minimal SU(5) model as we illustrate in the following section. 
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4. SU(5) example 



In minimal supersymmetric SU(5) [Tq] the GUT breaking is due to the VEV of an adjoint 
superfield E € 24. The superpotential for E includes the soft terms Ay, By as follows: 

= M s tr E 2 (l - B^6 2 ) + ^A E tr E 3 (l - A s 6 2 ) . (4.1) 

6 

The effect of Ay and By, is to give an F-term to (E) and to shift its magnitude by a 
small amount: 



(E) = v s 



1 + (A E - 5 S )0 2 



Ay , vs = 8\/l5 



As 



1 + 



As -^s 
2M S 



/ — My 
~8VT5-^, (4.2) 
As 



where Ay = diag(2, 2, 2, -3, -3) />/60. 

The fermion multiplets in SU(5) are 10 and 5, and proton decay can proceed via the 
four field operators involving the combinations 10-10-5-5 or 10-10-10-10. At GUT scale 
the standard D=6 operator has coupling constant 

A 6 (G) = gl/M 2 A , with M\ = 5^/12 , (4.3) 

where 55 is the SU(5) gauge coupling. Using eq. Q2.6p we find the coefficients of the new 
operators: 



£ = By - A 1 



L0 



105 



^1010 



-m 



^l 2 



2mi + mt + \By - A Y 



(4.4) 



where m|, m 2 and m s are the soft masses of the 5, 10 fermion multiplets and of E itself. 

In the previous section we found that the new operators for proton decay are relevant 
when the squark masses are small while £ or to are larger, ~ lOTeV. From (4.4) we see that 
this may be realized when the soft mass m s or the analytic soft terms By — Ay are large. 
A large m s is not appealing, since ms enters the RG running of the Higgs soft masses 
and would induce large values for these, spoiling the picture of electroweak breaking. The 
same holds for Ay, since it also enters in the running of m s and other soft parameters 
(see e.g. p0[ ), and a large Ay would indirectly cause color breaking minima. On the other 
hand we note that By does not enter the evolution of other quantities and may be sensibly 
large, without driving all the other soft parameters to large values as well. 

The fact that soft £?-terms do not enter in any beta function is a general statement valid 
in the MS scheme, that follows from SUSY and the fact that in this scheme no spurious 
scales are introduced. It turns out that the soft i?-terms like By, being of dimension 
one, never enter any RG equation, at all loops. Specifically, this can be verified in RG 
equations of soft masses, for A-terms and for i?-terms themselves; for example we have, 
for the one-loop running of Ay and By from Planck to GUT scale: 6 



d 



63 



16^ -Ay = ^AyXy + 3A H \ Z H - 30 g£ m~ g5 



dt 

16vr 2 -^5s 
dt 



20 
21 

.10 



AYX 2 Y + 2A H X 2 H -20glmg 5 . 



(4.5) 
(4.6) 



The constants Ah, Ah are denned below, (4.8). To compare with evolution of other quantities see 
e.g. @. 
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Figure 4: Effect of dimension 4 and 5 on gauge mediated proton decay as a function of the 
soft SUSY-breaking terms in the heavy sector. The dot represents our reference value for the 
supersymmetric gauge mediated proton lifetime, taken to be 10 36 y. The shaded regions show 
current and 10-years expected limits on the p — > 7r°e + partial lifetime. Dashed lines mark the limit 
where B% start to affect the higgs soft masses and other soft parameters. 



From these equations one can also see that if a large By is generated at the Planck scale, 
it will not be substantially affected running down to the GUT scale. 

Of course one can not hope to raise one soft parameter without consequences. Even if 
By, does not appear in RG equations, going from the MS to a physical scheme, it will enter 
in finite corrections. In particular a large By will generate corrections to soft quantities [19, 



20], raising them as if actually SUSY were broken at the By scale. As we discuss below, 
this effect is relevant only when By exceeds ~ 50 TeV, with some model dependence. Below 
this limit the only physical effect of a large By is in the soft £ and u coefficients, where it 
can directly dominate in the D=5 and D=4 operators and thus enhance or even suppress 
the proton decay rate. 

To give a concrete estimate in the SU(5) example, we assume large By and use the 
soft coefficients £ ~ By and uj ~ |-Bs| 2 in eq. (|3.8| ), to find how the proton lifetime for a 
gauge-mediated channel like p — ► 7r°e + is modified: 



— 1 —1 i 

t soft — t susy I 1 + io TeV + 0,1 



By 



10 TeV 



(4.7) 



An explicit plot of the effect of large By is shown in figure [|, where we assume a 
reference value of 10 36 y for the proton partial lifetime. 7 We see that the effect can be 
rather evident: for example for negative By the proton decay can be made absolutely 
unobservable, while for By positive one can enter in the region of sensitivity of the next 
ten years water-cerenkov detectors [^] . We conclude that large soft SUSY-breaking terms 
for the heavy fields may significantly affect the proton decay rate even in the gauge mediated 
channels. 



7 This reference value corresponds to a Grand Unification scale of 2-10 16 GeV, and we remind that t si ^ sy 
scales as the fourth power of Mg, which is model dependent. 
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Let us now address the fine tuning problems related to the stability of the hierarchy 
in presence of a large By- This point is related to the problem of doublet-triplet splitting, 
as can be seen in the SU(5) example. 

The superpotential involves also the higgses H, H (transforming in the 5, 5) 

W(S, H, H) = W(T,) + M H {\ - B H 6 2 )HH + A H (1 - A H 6 2 )HY,H . (4.8) 



where W(S) is given in (|4.l| ). After the SU(5) breaking H, H leave the light MSSM higgs 
doublets H u , H^, with their soft masses mj, m 2 ,, and we get also the effective [i and B^ 
terms 

- B^ 2 )H u H d (4.9) 



where 



H = M H - ~4^\hvy , (4.10) 



fiB^ = -1=\ h vy(Ay -By-A h + B h ) + 0(By - Ay) 2 . (4.11) 
VoO 



Therefore two fine-tuning conditions are needed to achieve the electroweak symmetry break- 
ing at the correct scale, one for \i and another for In fact the mass matrix of the higgs 
scalars is: 

( ^ + ^1 nBp \ 
y liBp H 2 +m 2 J 

and all entries should be of the order of the electroweak scale. 

The second fine tuning ( 4.1 1| ) can be avoided by assuming universality of A and B terms 



(4.12) 



separately, Ay = Ah and By = Bh, as noticed in |23j| , so that (Ay — By — Ah + Bh) = 
and B^ is of the order of soft susy-breaking scale. In the case of large By ~ 10 TeV however 
one still gets a B^ term that is too large, therefore the right pattern of electroweak breaking 
can be obtained only by tuning the two independent parameters, the super symmetric \i 
and soft B^. 

Of course, this minimal SU(5) model is not realistic, and one should not be surprised to 
find that fine tunings are required. In the next section we describe how in specific models 
fine tunings can be avoided and one can have large 5-terms without spoiling the hierarchy. 

Before moving to more realistic models, we point out that generically there are also 
finite corrections induced by S-terms. For example in SU(5) By induces a shift of the 
analytic soft term for the higgses B^, of their soft masses and also of the gaugino masses. 
These corrections are loop suppressed: 

^ - ^ 2 By , Sm~ g ~ j^By , (4.13) 



with some model dependent numerical factors [22]. Since Ah — 55 — 0.7, the loop sup- 
pression factor is ~ 1/100, and we conclude that these corrections can be ignored as far 
as By < 50 TeV. Beyond this limit the gaugino mass and the higgs mass terms would 
need some fine tuning, to avoid breaking SUSY at a high effective scale or having an 
unacceptably large higgs mass. 
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5. Realistic models 



In minimal SU(5) model the problem of doublet-triplet splitting has only a technical solu- 
tion: fine tuning of fj,, eq. ( 4.1G| ) that is stable against radiative corrections; the situation 



ft ~ 



is then worsened by the need of another fine tuning in the soft terms, eq. fl4.1l|) . 

Moreover, to achieve the right electroweak scale of order 100 GeV with a large B : 
lOTeV would require a fine tuning with the fi term, while there is no apriori correlation 
between these two parameters. 

This problem gets another twist in realistic models in which the doublet-triplet split- 
ting problem is solved without fine tuning. In particular, in SU(5) this can be done via the 
"Missing Doublet Mechanism" (MDM) Q, and in SO(10) via the "Missing VEV Mech- 
anism" (MDM) while in SU(6) via the "Goldstones instead of Fine Tuning" (GIFT) 
Mechanism |26[| . 

In particular, in all these models the soft parameters like By or Ay for the heavy GUT 
breaking superfields can be taken much larger than that of matter superfields, without 
creating additional fine-tuning problems. Let us briefly describe them here. 



In SU(5), the missing doublet model [24] contains the Higgs superfields in represen- 



tations ~ 75, H ~ 5, H ~ 5, \P ~ 50, \P ~ 50, with the following superpotential 
terms: 

W = M$ 2 + A$ 3 + Mi*§ + XiHm + \ 2 5m + fiHH (5.1) 

with M and Mi being the mass parameters order Mq and A's being the order 1 coupling 
constants. SU(5) is broken to SU(3)xSU(2)xU(l) by the VEV of <E> which also generates 
the mixing between the color triplet fragments in the Higgs 5- and 50-plets, whereas there 
are no doublets in the 50-plets. In this way, all color triplets are heavy, with mass of order 
Mq, while the doublets in H, H remain light, with mass given by the /i-term. Obviously, 
in this theory the soft parameter B§ can be taken large without inducing a large B^ (still 



inside the limits set by the induced finite corrections like ( 4.13 )). 



For SO(10), in the missing VEV model [25], the philosophy is similar: the Higgs dou- 



blets remain massless because the GUT-breaking fields have zero VEV along the direction 
that would give them a mass, whereas it couples to the triplets with non-zero VEV. There- 
fore also in this case the protection of the doublet sector is due to group theoretical reasons, 
therefore large soft terms ~ lOTeV in the heavy sector will not influence /i and B^ and 
the electroweak scale will not be destabilized. 



In the SU(6) model [26], the SU(6) gauge symmetry is broken by two sets of super- 
fields: one contains an adjoint representation S ~ 35, that leads to the breaking channel 
SU(6)^SU(4)x SU(2)x U(l), and the other contains two fundamental representations 
H ~ 6 and H ~ 6 that break SU(6)— >SU(5). As a result the two channels together 
break the SU(6) gauge symmetry down to SU(3) xSU(2) xU(l). Also, one assumes that 
the Higgs superpotential does not contain the mixed term HUH, so that it has the form 
W = W(E) + W(H, H), where 

W(E) = MY? + AS 3 , W(H, H) = Y(HH - V 2 ) . (5.2) 
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As a result there is an accidental global symmetry SU(6)exSU(6)h, which independently 
transform £ and H, H superfields. Then, in the limit of unbroken supersymmetry the 
MSSM Higgs doublet H u , appear as massless goldstone superfields built up as a combi- 
nation of doublet fragments from £ and H, H, that remain uneaten by the gauge bosons. 
Therefore in this limit \i vanishes exactly. 

Supersymmetry breaking terms like A-£, By, shift the VEVs and also give F-terms to 
them, therefore generating term for the MSSM Higgses. However, since these terms 
also respect the global symmetry SU(6)sxSU(6)h, the mass matrix of the Higgses ( 4.1 2| ) 
is degenerate and so one Higgs scalar (combination of the scalar components of H u and 
Hd) still remains massless. Thus, even with arbitrary B^ that give ~ B^ ~ B-£, there 
is an automatic relation between [i and B^ terms that guarantees that the determinant 
of ( 4.1 2| ) vanishes. 

This degeneracy is removed only by radiative corrections due to Yukawa terms that 
do not respect the global symmetry, and the resulting Higgs mass will be of the order 
of /j, and B^, given by the mismatch in their renormalization. Therefore, in the case of 
large B-£, ~ lOTeV we are still left with a "little" hierarchy problem of the electroweak 
scale stability against lOTeV. However by enlarging the gauge symmetry this issue can 
be avoided. In fact one can have that lOTeV is only an intermediate scale where an 
extra global symmetry guarantees the protection of the electroweak scale, the so called 
super- little-higgs mechanism [28|. 



6. Conclusions 

In this paper we have studied the effects of soft SUSY-breaking terms on proton decay in 
SUSY GUT theories. While the dominant effect in SUSY GUT comes from D=5 higgs- 
mediated operators, these are very model dependent and may be suppressed by specific 
constructions. Here we have focused on gauge mediated effective operators, that are usually 
unavoidable. 

We have shown how soft terms enter into the gauge-mediated effective operators for 
proton decay: while the supersymmetric operators are of dimension 6, SUSY breaking 
always induces new operators of dimension 5 and 4. 

We computed their renormalization from the GUT to the SUSY scale, that amounts 
to a small mixing of the D=6, 5, 4 operators through the gaugino masses. 

The new operators are dressed via gaugino exchange and transformed into D=6 four- 
fermion operators, and have the same suppression factor Mq 2 of the standard D=6 oper- 
ators. They however have numeric coefficients that depend on the ratio of soft-breaking 
parameters in the heavy and light sectors. 

When all the soft breaking parameters are of the same order, the dressing loop factors 
are small enough to suppress these new operators. However, we note that the 5-terms 
in the heavy-Higgs sector may be substantially higher than the standard soft masses, and 
they do not mix with soft masses under renormalization. Finite corrections are present 
which are irrelevant when the heavy 5-terms are smaller than ~ 50TeV. 
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The heavy higgses soft-terms then enter the GUT breaking process and lead to ob- 
servable effects on D=6 proton decay. £>-terms as low as lOTeV can lead to substantial 
effects on the proton decay and, depending on their sign, may enhance or even suppress 
the proton decay rate in gauge mediated channels. 
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